Spectral averaging techniques 
for Jacobi matrices with matrix entries 



Christian Sadel, Hermann Schulz-Baldes 

Department Mathematik, Universitat Erlangen-Niirnberg, Germany 

Abstract 

A Jacobi matrix with matrix entries is a self-adjoint block tridiagonal matrix with in- 
vertible blocks on the off-diagonals. Averaging over boundary conditions leads to explicit 
formulas for the averaged spectral measure which can potentially be useful for spectral anal- 
ysis. Furthermore another variant of spectral averaging over coupling constants for these 
operators is presented. 



1 Introduction 

Many variants of spectral averaging for one-dimensional Sturm-Liouville or Jacobi operators 
are known |CLj . If such operators depend on some continuous parameters, then the spectral 
averaging principle states that the spectral measures averaged over these parameters with respect 
to a measure with density are themselves absolutely continuous. In refinements useful for a 
detailed spectral analysis, it is possible to prove that they are even equivalent to the Lebesgue 
measure |dRT[ IdRMS] . The continuous parameters are typically boundary conditions or coupling 
constants. 

For Jacobi operators with matrix entries the only contribution seems to be due to Carmona 
and Lacroix [CLj . Unfortunately, their work does not give all the details of proof and the 
presentation is not conceptually structured nor does it cover full generality. Part of this work, 
in particular Theorem [2], is thought to fill these gaps. The main ingredient of the proof is the 
Cauchy formula for integration over the unitary group as proven by Hua pua ] (it is recalled 
in an appendix). Theorem [2] leads to a formula (stated in Theorem |3]) establishing a close link 
between spectral properties of the Jacobi operators in the limit point case and their formal 
solutions expressed in terms of the transfer matrices. We believe that Theorem [3] can potentially 
be a useful alternative tool (other than Kotani theory |KS] ) for proving existence of absolutely 
continuous spectrum. Finally Theorem H] provides a matrix version of a well-known identity of 
rank one perturbation theory, showing that averaging of the spectral measure over both boundary 
conditions leads to the Lebesgue measure. As an application, spectral stability results w.r.t. 
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local perturbations are presented. It is shown that also averages over fewer parameters than the 
whole set of boundary conditions lead to averaged spectral measures which are equivalent to the 
Lebesgue measure, at least locally in energy. This last part generahze the results in [dRMSj . As 
we lack a subordinacy theory for Jacobi matrices with matrix entries the applications to spectral 
theory of that paper do not carry over. 

This work complements our prior works [SBl] on Sturm-Liouville oscillation theory and |SB2j 
on Weyl theory for Jacobi matrices with matrix entries, so the basic notations and setup are 
chosen accordingly. Heavy use is being made of the matrix Mobius transformation on which 
there is an abundant literature (see the references in [SBH ISB2j ). but the main facts relevant 
for the present purposes are resembled in an appendix and all their short proofs are given in 
[SBTllsfe] . 

Acknowledgment: This work was supported by the DFG. We also thank the Newton Institute 
for hospitality. 

2 Setup and review of needed results 

Notations: The matrix entries of the Jacobi matrices are of size L G N. Matrices of size L x L 
are denoted by roman letters, those of size 2L x 2L by calligraphic ones. The upper half-plane 
Ui is the set of complex L x L matrices satisfying i{Z* — Z) > 0. Its closure U^, is given by 
matrices satisfying i{Z* — Z) > 0. The boundary is a stratified space dU^ = U^^-^diU^, where 
dilJi contains those matrices in Ul for which the kernel oi Z* — Z is /-dimensional. 

2.1 The Jacobi matrix and its resolvent 

Fix two integers L, N & N and let (Tn)n=2,...,N and {Vn)n=i,...,N be sequences of respectively 
invertible and self-adjoint L x L matrices with complex entries. Furthermore let the left and 
right boundary conditions Z and Z be also self-adjoint L x L matrices. Then the associated 
Jacobi matrix with matrix entries ^ is by definition the self-adjoint operator acting on states 

= (0n)„=i,...,ive£2(i,...,iV)®C^by 

{H^^^ (f))n = T„+10„+1 + Vn(f)n + T*(f)n-l , n = 1, . . . , , (1) 

where Ti = T^+i = 1, together with the boundary conditions 

00 = Z (f)i , = - Z (pN ■ (2) 

If Z = and Z = 0, one speaks of Dirichlet boundary conditions at the left an right boundary 
respectively. It will be useful to allow also non-selfadjoint boundary conditions Z, Z G Ul hence 
giving rise to a possibly non-selfadjoint operator . One can rewrite as an NL x NL 
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matrix with L x L block entries: 



H 



N 



( Vx-Z T2 

T* K 



v 



(3) 



At times, our interest will only be in the dependence of the right boundary condition Z, and 
then the index Z will be suppressed. 

As for a one- dimensional Jacobi matrix, it is useful to rewrite the eigenvalue equation 

{Hl^(t))n = zK. n = l,...,Ar, (4) 
for a complex energy z G C in terms of the 2L x 2L transfer matrices defined by 



n 



T- 







= 1, 



namely 



^n+l 



+lV^n+l 



T 



nYn 
'n-1 



n = l,...,N . 



(5) 



(6) 



This gives a solution of the eigenvalue equation (jl]) which, however, does not necessarily satisfy 
the boundary condition ([2]). Now z G C is an eigenvalue of ^ if and only if there is a solution 
of (jl]), that is produced by which satisfies ([2]). As is well-established, one can understand 
(12]) as requirement on the solution at sites 0, 1 and A^, + 1 respectively to lie in L-dimensional 
planes in C^'^. The corresponding two planes are described by the two 2L x L matrices (one 
thinks of the L columns as spanning the plane) 



(7) 



Then the boundary conditions ([2]) can be rewritten as 



4>o 



G <I^C^ 



Tn+i(Pn+i 

4>N 



G ^7X^ 



One way to attack the eigenvalue problem is to consider the L-dimensional plane $^ as the initial 
condition for an evolution of L-dimensional planes under the application of the transfer matrices: 



Z ■ 



(9) 



Because the transfer matrices are invertible, this produces an L-dimensional set of solutions of 
(in])- With the correspondence 



T 



n+l Vn+1 
0: 



(10) 



this also gives a matricial solution (;/>^ of (jlj). Due to the initial condition in (Q the left boundary 
condition at sites 0, 1 is automatically satisfied. The dimension of the intersection of the plane 
$^ with the plane $2 gives the number of linearly independent solutions of (jlj) at energy z, and 
therefore the multiplicity of z as eigenvalue of . 

Given ([9]), but also its own sake, it is natural to introduce the transfer matrices over several 
sites by 

T\n,m) = T„^ ■ . . . ■ T^+i , n > m , (11) 

as well as T^{n^ n) = 1 and T^{n, m) = T^{m, n)~^ for n < m. With this notation, the solution 
of the eigenvalue equation (jl]) satisfies $^ = T''{n,m)^^ and, in particular, $^ = T^(n, 0)$^- 
Of particular importance will be the transfer matrix T^{N,0) across the whole sample. Let us 
introduce the notations 

riN,0) = , (12) 

where all entries are L x L matrices. These matrices will intervene in many of the results below. 
Let us point out that T^{N,0) and all its entries do not depend on the boundary conditions Z 
and Z. The transfer matrix including boundary conditions is then 

1 z\ ( 1 o \ _ ( ai, + zc%-bIjZ-zdIjZ bIj + zdi,\ 

Now we introduce the resolvent. Let 7r„ : C'^ C^'^ for n = 1, . . . , N denote the partial 
isometry 

7ln\l) = \nj) , l = l,...,L , 

where the Dirac notation for localized states in ® is used. Then the L x L Green's matrix 
is given by 

G^(Z,Z,n,m) = <(i7f ^ - ^ 1)~V^ . 



Proposition 1 [gB2] For Z,Z e Vl, 



G%iZ,Z,l,l) 



+ ^ B%Z — ZD%Z 



1 -1 



[S^ + ZD 



+ ZCf^ — BIjZ — ZDfjZ 



'N' 



[Bl + ZDl] 



2.2 Parametrization of the boundary conditions 

The underlying hermitian symplectic structure is an important ingredient in most of the equations 
of the last sections, in particular in their proofs. It is necessary in order to understand what the 
adequate spectral averaging over the boundary conditions is. This section first recalls basic fact 
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about the symplectic structure, which will then be applied below. Let the symplectic form J be 
the 2L X 2L matrix defined by 

'0-1 
1 



J 



An L-dimensional plane described by a 2L x L matrix $ of maximal rank is called Lagrangian 
(or also isotropic, or simply symplectic) if = 0. 

Two L-dimensional planes described by 2L x L matrices $ and \l/ are called equivalent if 
there exists c G G1(L, C) with $ = \E'c. The Lagrangian Grassmannian is by definition the 
set of equivalence classes of Lagrangian planes. It is difficult to track back the original reference 
for the following result (it probably predates [Bot j ) . A short proof can be found in [SBl] where 
it is also shown how two natural symmetries are implemented. 

Proposition 2 The Lagrangian Grassmannian hi is identified with the unitary group U(L) via 
the real analytic diffeomorphism U : ^U(L) given by 



Um^) = {a-ib){a + zbr\ ^ = [l) 



Due to this theorem there is a natural measure on the Lagrangian Grassmannian given 
by the pull-back under n of the Haar measure on the unitary group. 

The Lie group conserving the (hermitian) symplectic structure is the (hermitian) symplectic 
group SP(2L,C) defined by those 2L x 2L matrices T satisfying T*JT = J. Clearly, if $ 
describes a Lagrangian plane, then so does T$ for any T gSP(2L, C). Isomorphic to the 
hermitian symplectic group is the Lorentz group of signature (L, V) defined by U(L, L, C) = 
CSP(2L, C)C*, where C is the Cayley transformation introduced as the matrix 

c 1 ^1 -'1 



72 V 1 '1 

Next let us exhibit explicitly the symplectic structure in the equations of Section 12.11 Both 
of the planes $^ and $^ used as boundary conditions in (IHl) are Lagrangian in the above sense. 
Actually there are many Lagrangian planes which cannot be written in this way, but they form 
a set of zero measure. Due to Proposition [2l it is natural to identify the left and right boundary 
conditions with unitary matrices: 

ij = n([|.^].) , u = uii^zU . (14) 

In other terms, this means U = C ■ {—Z) and U = —C ■ Z. Furthermore, let us set 

= n([$f].) , 

where is the solution which automatically satisfies the left boundary condition. Then 
Uq = tj . Using the correspondence (IT^ . we also set 

Gl,{U,U) = G^(Z,Z,1,1). 
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As z ^ Gff{U, U) G Ul is analytic in z for '^'m{z) > 0, the Herglotz representation theorem 
|GT] associates a matrix-valued (averaged spectral) measure: 



2.3 The oscillation theorem 

The oscillation theorem is another application of the parametrization of boundary conditions. It 
is stated for sake of completeness and because it will be used in the proof of the result of Section HI 
Due to Section E^l £^ e R ^ $f is a path of Lagrangian planes and for each E the dimension 
of its intersection with the right boundary condition is the multiplicity of E as an eigenvalue 
of ^. This intersection number was introduced by Bott [Botj precisely for the study of the 
eigenvalue calculation of Sturm-Liouville operators, the continuous analogues of Jacobi matrices. 
Later on it was rediscovered by Maslov and a detailed survey of its properties is included in [SBl]. 
The intersection number can be conveniently calculated using the associated unitary and this 
leads to the following theorem which was proven in |SBlj under the supplementary hypothesis 
that the T„'s are positive, but the proof directly transposes to the slightly generalized situation 
considered here. 

Theorem 1 Let E ^ M., N > 2, and {for sake of simplicity) the right boundary condition be 
Dirichlet, that is Z = 0. Then there are L strictly increasing real analytic functions 6§ i : ^ M., 

I = 1,...,L, such that e*^^-' are the eigenvalues ofU§. The multiplicity of E as an eigenvalue 
of is equal to the multiplicity of —1 as an eigenvalue of U§. Furthermore, the matrix 

7 {U§)*dEU§ is positive definite. 

2.4 Limit point operators 

If in the prior sections N = oo, then the right boundary condition Z is pushed to infinity. If 
this gives a well-defined (essentially self-adjoint) operator H^, one speaks of the limit point case. 
Various criteria for this can be given, the simplest one being that \\Tn\\ is uniformly bounded 
from below. Otherwise one needs the infinite operator having non-vanishing deficiency spaces 
and one has to consider various self-adjoint extensions. Here we restrict ourself to limit point 
operators. For these operators, limits 

G^iU) = hm GUU,U) , 

exist, are independent of U and are the Green function of H^- Its spectral measure is denoted 
by and obtained as the weak limit of p^,.- 
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3 Average over boundary conditions 



Let us write dU for the normalized Haar measure on U(L). 
Theorem 2 For '<sm{z) > 0, one has 

jdUG%{l,U) = [A% + iC'^]-'[B'^ + tDl,] (15) 
= [{Bir + ^{Dly] [{Air + z{Cir]-\ (16) 
Moreover, for E = ^e{z), and setting Qm{A) = ^(^ — A*) for any square matrix A, 

Jim [dUQm{Gl,il,U)) = [{Af,rAf, + iC^yC^]-\ (17) 

Remark 1 A formula similar to (ITTj) can be found in j CL] , but the latter authors use the average 
over the Haar measure on the symmetric space of symmetric unitaries instead of the group U(L) 
(moreover, their proof seems to have several gaps). 

Remark 2 It is easy to incorporate the left boundary condition U 1 using (fT3ll . 

Remark 3 One way to define the closed Weyl disc W|^r is as the image of the map Z E Ul ^ 
Gj^{l, Z,l,l). As proven in [SB2j the points in the Weyl disc can also be parametrized by 
6*^(1, Z, 1,1) = + (i?|^) 2 VF(i?^) 2 where and R]^ > are properly defined center and 
radial operators, and W e U(L) depends on Z, cf. |SB2j . Taking the average over W w.r.t. to 
the Haar measure in this representation immediately gives / dW G'}^^!, Z, 1,1) = S^, which is 
not equal to the r.h.s. of (ITS!) . The Jacobian of the change of variables Z does not seem 

to be known (nor be of great importance). 

Proof of Theorem [2l First let us note that one can use the Mobius transformation to express 
V = H([$^]^) = C ■ {— Z). Hence also Z = —C* ■ U. Starting from Proposition [H one therefore 
has 

j dUG'^{l,U) = jdU {A%-C*-UG;,r'{B'^-C*-UD%) . 

By 

j dUGlj{l,U) = jdU {Al,+ C*-UG'^)-'{B'^ + C*-UDl,) . 

(Alternatively to this argument, one could have defined the average on the l.h.s. by the r.h.s..) 
In order to be able to apply the Cauchy formula (12^ for Z = 0, it is sufficient to show the 
analyticity of the function 

f{Z) = {A% + C-ZClj)-\B% + C*-ZDl,), 

on the unit disc as well as its continuity on the closure D*^ (strictly speaking, one should 
consider the entries of the matrix- valued function /). This follows from Weyl theory jSB2j 
combined with the fact that (— C* ■ Z) is in the closed lower half plane (for Z G D^). The Cauchy 
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formula (!24|) for Z = now concludes the proof of (ITS!) because C* ■ = tl. Formula (fT6|) is 
proven similarly from the second identity in Proposition [1] 

It follows from the results of [SB2] that + zC^ = (A^(C^)-^ + 2l)C^ is invertible. 
Inserting 1 = [(A^ + 2C^)*]-i(^^ + ^C*^)* in dlSD shows 



As the transfer matrices at real energies are symplectic, the limit of vanishing imaginary part in 
the energy can be taken in this equation and that directly implies 

^lim [dUGUl,U) = [{A^^rA§ + {C§yC§]-' [{Af,rB§ + {C§yD§ + tl] , 

Qm(2)J,0 J 

where we used the identity {Aff)*D§ + {C§)*B§ = 1 holding for any symplectic matrix. The 
same calculation can be carried out starting from ( fT6i) and adding the results up gives ( fTTl) . □ 

One corollary of Theorem [2] is the following formula which links the averaged spectral measure 
defined by 

PO = J dUpl^ (18) 
to properties of the eigenf unctions of the transfer matrices at real energies. 
Corollary 1 For any Eq < Ei , 



\[p^^{[E,,E,])+p^^{{E,,E,))] = r^dE [|A^-<C*-(-?7)p+|C^-D^r-(-[/)p 

In the limit point case, the averaging in f|T8|) becomes irrelevant because ^ converges weakly 
to the spectral measure pjy of as — oo. This leads to the following formula for the spectral 
measure, which was already obtained by Carmona in the strictly one-dimensional case L = 1, 
and by Pearson |Pea] for one-dimensional Schrodinger operators. For sake of simplicity, let us 
set [> = 1 so that C* ■ {-&) = 0. 

Theorem 3 Let the semi-infinite Jacobi matrix H be in the limit point case. Then, for any 

Eq < El, 



^[pi([Eo,£;i])+pi((i?o,i?i))] = lim pdE [(A^)Mf + . (19) 



Next we also average over the left boundary condition If. The associated averaged spectral 
measure is equal to the Lebesgue measure, a fact also known from the case L = 1. 

Theorem 4 For any N and , 

4 j dU {U - 1)-^ P^{dE) (U* - = IdE . (20) 
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The proof of the theorem is based on the following integral identity. 
Lemma 1 Let < T e U(L, L) and V E U(L). Then 



dU 



-i -1 



Proof. Let / denote the integral appearing in the lemma. First let us use that < T G U(L, L) 

f W \ 

can be transformed into a normal form hy M = [ ^ ] e U(L, L) fl U(2L) = U(L) © U(L) 



where W, W G U(L), namely 



M*rM 



cosh(77) sinh(77) 
sinh(?7) cosh(?7) 



where r] = diag(r7i, . . . ,77^) is a diagonal matrix with non-negative entries. Let us denote the 
r.h.s. by 7^. Replacing this identity, one obtains 



{w'vy / du 



wu(w'vy 
1 



wuiw'vy 



-1 



(WV) , 



so that using the invariance of the Haar measure dU one realizes that it is sufficient to consider 
the case V = 1 and T = 7^. First suppose 77 > 0. In this case, sinh(77) is invertible and 



dU [cosh.{ri)U + U cosh(?7) + U smh.{'r])U + sinhr/] U 



dU 



U 



cosh(77) + 1 
sinh(?7) 



-1 



(sinh(77))"^ ( U + 



cosh(?7) — 1^ ^ 
sinh(?7) 



U 



To simplify notations, define a = ^"g.^^^^^"^"*" and (3 = Note that < /3 < 1 < a, a[3 = 1 



dU{U + a)-^U = 0, 



and {a — (3) sinh(r7) = 2. Now using 

f dU{U + l3)-^U = I dU^{-U-^f3Y 

n>0 

one gets with the resolvent identity 

/-I = 2^ (i[/[((?7 + a)-^-(«-/9)-^)sinh(r7)-i([/ + /3)-^[/] = - ^ dU{U + a)-^U = 0, 

and the proof is complete in the case i] > 0. By continuity of the integral one also recovers the 
case 1] > 0. □ 

Proof of Theorem HI This is based on Corollary [H Actually it is sufficient to show that the 
average of the integrand on the r.h.s. of Corollary [1] satisfies 

4 [dU{U-l)-^ \\Af,- B§C*-{~U)\^ + \C^- D§C*-{-U)\^Y\u* = 1. 



(Note that the various inverse appearing in this formula do not exist on a set of zero measure.) 
This identity reduces to 

-1 



dU 



U 
1 



CT^{N, 0)*r^(A^,0) c* 



u 
1 



1 . 



Now r^(A^,0)*T^(A^,0) is a positive symplectic matrix and hence CT^(iV, 0)*T^(A^, 0)C* is a 
positive matrix in U(L, L). Hence Lemma [1] shows that this identity indeed holds. □ 



4 Spectral measures averaged over coupling constants 

As an application of the results of Theorem [3] we consider here a particular one-parameter family 
of Jacobi matrices with matrix entries obtained by a local positive perturbation and show that 
the associated averaged spectral measure is under certain conditions absolutely continuous. 

Let if be a Jacobi matrix with matrix entries in the limit point case, and let tt^ : (C^)^ 
denote the partial isometry onto the nih site. For real positive semi-definite matrices {Wn)i<n<N 
and G M define 

N 
n=l 

Furthermore, let [jj) be the finite Jacobi matrix obtained by projecting H[n) to the first 
sites. It is of the form ([3]) with Vn replaced hj + pWn and Z = Z = 0. Because the 
perturbation H{ij,) — H is increasing in /i, the eigenvalues of H{fj,) are increasing functions of 
Finally let be the matrix-valued spectral measure of H{ij,) and define the averaged spectral 
measure corresponding to an interval / = [fiQ, ni] by 

p = j^dn Tr(p(^)) . 

Theorem 5 Suppose that Wn > and Wn+i > for some n = 1, . . . , N — 1. Let I = [jj,Q, ^i] 
be sufficiently large such that there are 2L eigenvalues of H^{ii) passing by E as fi varies in I. 
Then p is equivalent to the Lebesgue measure in a neighborhood of E. 

As can be seen from the proof below the hypothesis can be somewhat relaxed. For L = 1 
the result was proven in |dRMS] . Similar as in |dRMS] . the condition on the size of I can 
also expressed in terms of an associated Birman-Schwinger operator and, furthermore, it is also 
possible to consider several parameter spectral averaging instead of over just one parameter ^. 
On the other hand, the applications to spectral analysis do not carry over immediately, because 
the subordinacy theory is not yet developed for Jacobi matrices with matrix entries. 

First we need to fix some notations. Just as H{^), all objects of the previous sections 
depend on a supplementary parameter /i. In particular, we will write T^{^) and T^{n,m, 
Furthermore, let us introduce the Dirichlet solutions = T^(A^, 0, /i) (J) and the matrix 

Pn{p) = -n'^ipTJ^^Tip)- 

The proof of Theorem O will be based on the following criterion. 
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Lemma 2 Suppose: 

(i) There exist positive constants Ci, C2 such that Ci 1 < P§{li) < C2 1 for all & I. 
(ii) 

l^^rnd, log(n«^(/.))) < -L. 
Then]) is equivalent to the Lebesgue measure in a neighborhood of E. 

Proof. We write p{Eq,Ei) for \ [p{[Eo,Ei]) + -p{{EQ, Ei))]. Let us start by integrating ffTg]) 
over ^ and using the dominated convergence theorem as well as Fubini's theorem: 



p{Eq,Ei) = lim 



Eo 



dE / dp Tr 



T^(M,iV)vl/;^^(/x) 



Now for positive semi-definite operators A, B with < Cil < B < C2I, 

Ti (AB) < Tt{A) < Tt{AB) . 
Applying these bounds for B = P§{p) shows that 



lim 



dE j dp Tr 



r%M,N)^-f{p) p^{p) 



where the sign ^ means that we have two-sided bounds. As T^{M, NyjT^{M, N) = J, 
^^{p) = T^{M, A^)\l/^^(/i) and T^{M, N) does not depend on p, this can be rewritten as 



— lim 



Eo 



dE / dp Tr 



— TT lim 

Af^oo 



El 



Eo 



dE l^^^rnd, log(n(vl/-f (^))) 



where the second identity is checked in |SBlt Lemma 4]. Now the expression under the integral 
J dE on the r.h.s. is precisely the pairing J-^ut of the Arnold cocycle to with the path r(^) = 

^m^(a*) ~ ^'^(^) ^)^)v^(/^)' P ^ I, the Lagrangian Grassmannian (actually here this is 
a path in the real Lagrangian Grassmannian because H{p) is real). Hypothesis states something 
about the pairing with the path T'{p) = namely J^, lo < —L. However, these two 

paths are related by the multiplication with the symplectic matrix T^{M, N). Hence by |SBH 
Proposition 4]: 



I- 




< L 




Jr 





Therefore < C3 < — Jp < C4 where the upper bound follows from compactness of / and the 
constants are independent of M. Replacing this shows 



p{Eo,Ei 



lim 

M^oo 



El 



dE 



El — Eq 



Eo 
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which is precisely the claimed equivalence of p with the Lebesgue measure. □ 
Proof of Theorem [5l First of all, 

N 

T'^{N,0,^rj*d,r'^{N,0,^) = J2 {n- 1,0, fiYT^^if,)* J* {d,T^^ if,)) T^{n- 1,0, fi). 

n=l 

But 

r^'^if^yrd^r^'^ip) = 



(T-i)W„T-i 




is positive semi-definite, and the arguments in the proof of [SBlt Proposition 6] show that 

( (^»«)"'^^-^- I ) + TIM- ( J ^ „ 

(strict positivity), whenever Wn > and Wn+i > 0. In the latter case the above sum is therefore 
strictly positive. Hence the hypothesis of Theorem imply P§{id) > Ci > 0. Compactness of 
/ thus imply that hypothesis (i) of Lemma [2] holds. Hypothesis (ii) follows from the oscillation 
theorem as stated in Theorem[TJ In fact, as fi increases each phase 0§i{fi) decreases. The integral 
in (ii) is the total phase (sum of all 6'^;(/i)'s, in units of 2tt) accumulated as fi varies in /. If 
L + K eigenvalues pass by as /i varies, the total phase has to change by at least K. Hence the 
hypothesis of Theorem [5] imply (ii) of Lemma [2l □ 



Appendix A: reminder on Mobius transformations 

This appendix resembles the basic properties of the Mobius transformation as they are used in 
the main text. All proofs are contained in |SBlj . Complex matrices of size 2L x 2L are denoted 
by mathcal symbols, those of size L x L by roman letters. 

The upper half-plane and unit disc (also called Cartan's first classical domain) are defined by 

Vl = {Z e Mat(L X L, C) | i{Z* - Z) > } , Bl = {U e Mat(L x L,C) \ U*U < 1 } , 

where y > means that Y is positive definite. If Z E Vl, then Z is invertible and —Z^^ G Vl- 
Moreover, for any V = V* and any invertible T, one has Z + V G Vl and T*ZT G U^. The 
formulas 

U = {Z - il){Z + il)-^ , Z = i{l + U){1 - U)-^ , (21) 

establish an analytic diffeomorphism from Ul onto Dl- The boundary 83 l of Dj;, is a stratified 
space given as the union of strata S^D^, I = 1, . . . , L, where di^^ is the set of matrices U for which 
U*U < 1 and rank(l — U*U) = L — I. By Proposition [2] the maximal boundary 9^©^ = U(L) 
is identified with the Lagrangian Grassmannian h^. Similarly, the boundary of is stratified, 
but this will play no role here. 

The Mobius transformation (also called canonical transformation or fractional transformation) 
is defined by 

r-Z = {AZ + B){CZ + D)-\ r = ( ^ ^ ) G GL(2L,C), Z G Mat(LxL,C) , (22) 
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whenever the appearing inverse exists. For T as in (122|) and as long as the appearing inverse 
exists, the inverse Mobius transformation is defined by 

W -T = [WC - A)-'^ {B - WD) , W G Mat(L x L, C) . (23) 

The Mobius transformation is a left action, namely {TT') ■ Z = T ■ (T' ■ Z) as long as all objects 
are well-defined. It is well-known that if T gSP(2L,C) and Z G Ul, the Mobius transformation 
T ■ Z is well-defined. 



Appendix B: Cauchy formula for Cartan's classical domain 

The results of this section are proven in |Hua] . Let 3^ be the topological closure of 3^ and let 
dU be the normalized Haar measure on its maximal boundary = U(L). 

Theorem 6 For any continuous function f : 3^ C which is analytic on 3^, one has for all 

z e3L, 

f{Z) = [dU det(l - ZU*)-"- fiU) ■ (24) 



For the proof of Theorem [2] we only need the case Z = 0. Moreover, only intervene functions 
/ which are of the form f{U) = F{C* ■ U) where C* -U is the Cayley transform of a unitary matrix 
(written in the notations of the appendix) and hence hermitian, and F is a complex function on 
the hermitian matrices. The Cayley transform C* ■ U is not defined for all unitaries U, but it is 
defined on a set of full measure. The change of variables formula to the Lebesgue measure d^ on 
(real and imaginary part of each entry of) the Hermitian matrices is now 

I dU F{C -U) = c j di det(l + CiY'^ F{i) , 

where c is a normalization constant (which is given in |Huaj ) . As the measure d^ det(l + ^*^)~^ 
is invariant under the refiection ^ — ^, it follows that 

j dU F{C* -U) = J dU F{- C*-U) . (25) 
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